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Abstract 



It is proved that the system of string equations of the dispersionless 2-Toda hierarchy which 
arises in the planar limit of the hermitian matrix model also underlies certain processes in Hele- 
Shaw flows. 



1 Introduction 

> 

00 

I matrix models [I] -[3]. Thus, the partition functions 

d 



o 



The Toda hierarchy represents a relevant integrable structure which emerges in several random 

s, the partition functions 

Zjv (Hermitian) = J dH exp (tr(^t fc tf fc )), (1) 



k>l 



Z N (Normal) = ! dM dM f exp (tr (M M f + ^(t k M k + t k M f fc ))) , (2) 

of the hermitian (H = H^) and the normal matrix models ([M, M*\ = 0) , where N is the matrix 
dimension, are tau-functions of the 1-Toda and 2- Toda hierarchy, respectively. As a consequence 
of this connection new facets of the Toda hierarchy have been discovered. Thus the analysis of the 
large A-limit of the Hermitian matrix model lead to introduce an interpolated continuous version of 
the 2-Toda hierarchy: the dispersionful 2-Toda hierarchy (see for instance [4]). On the other hand, 
the leading contribution to the large A-limit (planar contribution) motivated the introduction of a 
classical version of the Toda hierarchy [4] which is known as the dispersionless 2-Toda (d2-Toda) 
hierarchy. 

Laplacian growth processes describe evolutions of two-dimensional domains driven by har- 
monic fields. It was shown in [5] that the d2-Toda is a relevant integrable structure in Lapla- 
cian growth problems and conformal maps dynamics. For example, if a given analytic curve 
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7 (z = z(p), \p\ = 1) is the boundary of a simply-connected bounded domain, then 7 evolves with 
respect to its harmonic moments according to a solution of the d2-Toda hierarchy. These solutions 
are characterized by the string equations 

z = m, m = —z. (3) 

Here (z, m) and (z, fh) denote the two pairs of Lax-Orlov operators of the d2-Toda hierarchy. As 
it was noticed in [5] -[9], this integrable structure also emerges in the planar limit of the normal 
matrix model (|2]) and describes the evolution of the support of eigenvalues under a change of the 
parameters tk of the potential. 

The present paper is motivated by the recent discovery [10] of an integrable structure provided 
by the dispersionless AKNS hierarchy which describes the bubble break-off in Hele-Shaw flows. In 
this work we prove that this integrable structure is also characterized by the solution of a pair of 
string equations 

z = z, m = m, (4) 

of the d2-Toda hierarchy. Since the system describes the planar limit of ([1]), it constitutes a 
common integrable structure arising in the Hermitian matrix model and the theory of Hele-Shaw 
flows. 

Our strategy is inspired by previous results |llj-[12j on solution methods for dispersionless 
string equations. We also develop some useful standard technology of the theory of Lax equations 
in the context of the d2-Toda hierarchy. 

The paper is organized as follows: 

In the next section the basic theory of the d2-Toda hierarchy, the method of string equations 
and the solution of (JH) are discussed. In Section 3 we show how the solution of (Jl|) appears in the 
planar limit of the Hermitian matrix model and the Hele-Shaw bubble break-off processes studied 
in P]. 



2 The dispersionless Toda hierarchy 
2.1 String equations in the d2-Toda hierarchy 

The dispersionless d2-Toda hierarchy[J] can be formulated in terms of two pairs (z,m) and (z,rh) 
of Lax-Orlov functions, where z and z are series in a complex variable p of the form 

u\ v 

z = p + u-\ 1 , z = - + v + vip-\ , (5) 

p p 

while m and fh are series in z and z of the form 

3=1 3>1 3=1 3>1 

The coefficients in the expansions ([5]) and (|6]) depend on a complex variable x and two infinite sets 
of complex variables t := (ti, t%, . . .) and t := (tx, . . .). The d2-Toda hierarchy is encoded in the 
equation 

00 

dz A dm = dz A dm = d^logpdx + ^ {{z j )+ dtj + diA^j. (7) 

3=1 



2 



Here the (±) parts of p-series denote the truncations in the positive and strictly negative power 
terms, respectively. As a consequence there exist two action functions S and S verifying 

oo 

dS = mdz + logpdx + ^ {[z j ) + dtj + dt^j, 



dS = mdz + logpdx + ^ l(z j ) + dtj + (z j )_ dt 



i=i 

and such that they admit expansions 

3=1 J ■ J=l 3>1 

From ([7]) one derives the d2-Toda hierarchy in Lax form 

where JC = z, m, z, fn, and we are using the Poisson bracket {/, g} := p(f p g x — fxdp)- 

The following result was proved by Takasaki and Takebe (see [1]): 

Theorem 1. Let (P(z,m),Q(z,m)) and (P(z,m),Q(z,m))) be functions such that 

{P,Q} = {z,m}, {P,Q} = {z,m}. 

If (z,m) and (z,m) are functions which can be expanded in the form ©-(H]) and satisfy the pair 
of constraints 

P(z,m) = P(z,m), Q(z,m) = Q(z,m), (10) 

then they verify {z, to} = {z,m} = 1 and are solutions of the Lax equations ([9]) for the d2-Toda 
hierarchy . 



Constraints of the form (I10p are called dispersionless string equations. In this paper we are 
concerned with the system @ . The first equation z = z of defines the 1-Toda reduction of the 
d2-Toda hierarchy 

v . . 

z = z = p + u + -, (11) 

p 

where 

u = d x S 2 , log v = -d x S - (12) 

As a consequence the Lax equations ([9]) imply that u and v depend on (t,t) through the combi- 
nation t — t. 

Due to (jlip there are two branches of p as a function of z 



p(z) = - ((z - u) + \J{z - u) 2 - 4vj = z 

(13) 



v 

u-- + 

z 



P{z) = -{{z-u)- y/(z- u) 2 - 4vj 



V 

- + 

z 



3 



To characterize the members of the dl-Toda hierarchy of integrable systems as well as to solve the 
string equations Q it is required to determine (z 3 )_(p{z)) and (z J ) + (p(z)) in terms of (u,v). By 
using (fTTj) it is clear that there are functions oy, fy, otj, which depend polynomially, in z such 
that 

d tj S(z) = (zi) + (p(z)) = oj + P jP (z), d- tj S{z) = (zJ)-(p(z)) = aj + Pjp(z), 

d t .S(z) = (Zi) + (p(z)) = Oj + (3jP(z), drS( Z ) = (zj)-(p(z)) = OCj + Pjp(z), 



and 



a j = z J -a j , Pj = -Pj- (14) 



Now we have 

ay + (3j p(z) = d tj S(z) = z> + O (~) , ay + f3 3 p{z) = d tj S(z) = -d tj S + O Q , (15) 
so that 

* j = l(zi-d tj S -(p + p)P J ), A"=(^) e . (16) 

where ( )® and ( ) e stand for the projection of z-series on the positive and strictly negative powers, 
respectively. Thus, by introducing the generating function 

= = 77 N9 i = k — ' r o = 1- (17) 

P-P ^{z - uf - 4v z k 



we deduce 



(**)+(p(z)) = z J -\ d tj S - ^ (zJ- 1 R) Q 



l(d tj S + r i) - J~ z (r j+ i ~ urj) +0[±). (18) 



Hence 



1 



^ S = -rj, d tj S 2 = - (r j+1 



ur 



3)i 



so that the equations of the dl-Toda hierarchy are given by 



9 tj u = - d x (r j+ i - u Tj) , d tj v = vd x rj. (19) 



Furthermore, we have found 

1 

Hence, the first terms of their asymptotic expansions as z — > oo are 



(**)_(p(z)) = -l rj + ^ (z^R) , = Ti + (**)-(?(*))• (20) 



(*>')-(?(*)) = J- (r i+1 -ur,-) + O(l) , (*>')+(p(*)) = ry + ^ (r i+a - ury) + (21) 
Notice that since ro = 1 and r± = u, these last equations hold for j > 0. 
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2.2 Hodograph solutions of the 1-dToda hierarchy 

In the above paragraph we have used the first string equation of Let us now deal with the 
second one. To this end we set 

CO CO 

m = m = Y^Jtj + J2^3 

3=1 3=1 

which leads to the following expressions for the Orlov functions (m,m) 



in 



(z) = Ytitj J- 1 + E-? (*i - *i) (^ _1 )-(pW). 

3=1 3=1 



(22) 



m(z) = Y^ih z j ~ x - Y.3 & - **) 

j=i i=i 



In order to apply Theorem 1 we have to determine u and v and ensure that (m, m) verify the 
correct asymptotic form (|5])-([6j). Both things can be achieved by reducing (|22p to the form 



x 

z 



1 CO 

+ E 4% = I>" & - ^) (^-(pC*)). 



i>2 



j'=i 



(23) 



1 CO 

7 + E ^ = -X>' ^ - tf-^+m), 



Z A — ' Z3 

i>2 



3=1 



and equating coefficients of powers of z. Indeed, from (|21f) we see that identifying the coefficients 
of z^ 1 in both sides of the two equations of ([23j) yields the same relation. This equation together 
with the one supplied by identifying the coefficients of the constant terms in the second equation 
of (|23p provides the following system of hodograph-type equations to determine (u, v) 



It can be rewritten as 



2 E^iife -^ 



o. 



(24) 



7 27ri ^(z - u) 2 - 4v 

dz z V 7 



0, 



(25) 



^ J 7 27ri y^z - u )2 - 4 V 



-2x, 



where 7 is a large enough positively oriented closed path and V z denotes the derivative with respect 
to z of the function 



V(z,t-t) :=E(*i-*J 
3=1 



))zK 



(26) 
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The remaining equations arising from (|23p characterize the functions Sj and Sj for j > 1 in 
terms of (u,v). Therefore we have characterized a solution (z,m) and (z,m) of the system of 
string equations verifying the conditions of Theorem 1 and, consequently, it solves the dl-Toda 
hierarchy. 

3 Planar limit of the Hermitian matrix model and 
bubble break-off in Hele-Shaw flows 

3.1 The Hermitian matrix model 

If we write the partition function (pQ) of the Hermitian matrix model in terms of eigenvalues and 
slow variables t := et, where e = 1/JV, we get 



„ n 

Z n (Nt)= H(dx k e NV ^))(A(x 1 ,--- ,x n )) 2 , V(z,t):=J2t k . 



z k . (27) 

The large iV-limit of the model is determined by the asymptotic expansion of Z n (N t) for n = N 
as N ^ co 

r N 

Z N (Nt)= / l[(dx k e Nv ^)\(A(x u --- ,x N )f, (28) 
JrN k=i 

It is well-known [3] that Z n (t) is a r- function of the semi-infinite 1-Toda hierarchy , then there 
exists a r-function r(e, x, t) of the dispersionful 1-Toda hierarchy verifying 

T(e,en,t)=Z n (Nt), (29) 

and consequently 

T (e,l,t) = Z N (Nt). (30) 
Hence the large iV-limit expansion of the partition function 

Z^(iVt)=exp(iV 2 F), F = ^J^F( 2fc ), (31) 

k>0 

is determined by a solution of the dispersionful 1-Toda hierarchy at x = 1. 

As a consequence of the above analysis one concludes that the leading term (planar limit) 
is determined by a solution of the 1-dToda hierarchy at x = 1. Furthermore, the leading terms of 
the ^-expansions of the main objects of the hermitian matrix model can be expressed in terms of 
quantities of the 1-dToda hierarchy. For example, in the one-cut case , the density of eigenvalues 

p{z) = M(z) ^J{z-a){z-b), 

is supported on a single interval [a,b]. These objects are related to the leading term of the 
one-point correlator [13] 

W(z) - 1 V ^(trMi) = - + —Y J- d^Z N {Nt) 
w W - N 2^ z i+x V rM > z + N 2 zj+ i dt 

j>0 3>l 3 
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in the form 

W i0) = ~V z (z)+mp(z). 
On the other hand, it can be proved (see for instance [14j ) that 



=m(z,l,t)- y £ i jt j zi- 1 , (32) 



so that ([20D and ([22]) yield 

1 oo 



i=i 



(33) 



OO ^ OO oo 

2 E 'i r ^-i - 2 E *i zj ~ l + 2 (P " P) E J ^ ( Z '" 2 * 

j=l j=l j=2 



Since we are setting tj = 0, Vj > 1, according to the first hodograph equation ([24"|) the first 
term in the last equation vanishes. Therefore the density of eigenvalues and its support [a, b] are 
characterized by 



p( z ):=±-( \ ) ^ {z - a){z - b) , 

2m \y/(z - a)(z - 6)/ e 

(34) 



a := u — 2 \/v, b := u + 2 yu, 
where we set x = 1 in all the x-dependent functions. Observe that according to ([33 



iirp(z) = ^V z {z) + - z +0[^), z^oo, (35) 
so that the constraint x = 1 means that the density of eigenvalues is normalized on its support 



\ p{z) dz = 1. 

J a 



Moreover, from (|25p we obtain 



dz Vz n I dz zV * o /o^ 

°> f -HZk -h w =Tf = " 2 > ( 36 ) 



7 2vri ^(^-0)^-6) 7 7 2vri ^ _ a )( z _ &) 

with 7 being a positively oriented closed path encircling the interval [a, b}. These are the qua- 
tions which determine the zero-genus contribution or planar limit to the partition function of the 
hermitian model 114 1-1 181. 
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3.2 Bubble break-off in Hele-Shaw flows 



A Hele-Shaw cell is a narrow gap between two plates filled with two fluids: say oil surrounding 
one or several bubbles of air. Let D denote the domain in the complex plane C of the variable A 
occupied by the air bubbles. By assuming that D is an algebraic domain [10], the boundary 7 of 
D is characterized by a Schwarz function § = S(A) such that 

A* = S(A), A G 7. (37) 

The geometry of the domain C — D is completely encoded in § and it can be conveniently described 
in terms of the Schottky double [9]: a Riemann surface 1Z resulting from gluing two copies H± of 
C — D trough 7, adding two points at infinity (00,00) and defining the complex coordinates 




In particular § dX can be extended to a unique meromorphic differential u on 1Z. 

The evolution of 7 is governed by D'Arcy law: the velocity in the oil domain is proportional to 
the gradient of the pressure. In the absence of surface tension, pressure is continuous across 7 and 
then if the bubbles are assumed to be kept at zero pressure, we are lead to the Dirichlet boundary 
problem 

'AV = 0, 011C-A 

V = on 7, (38) 

V ^ — log|z|, z — > 00. 

If one assumes D'Arcy law in the form v = —2 VP, then by introducing the function 

*(A):=f(A) + iP(A), (39) 

where £ and V are the stream function and the pressure, respectively, D'Arcy law can be rewritten 
as 

d t § = 2id x $, (40) 

where t stands for the time variable. 

In the set-up considered in [10] air is drawn out from two fixed points of a simply-connected air 
bubble making the bubble breaks into two emergent bubbles with highly curved tips. Before the 
break-off the interface oil-air remains free of cusp-like singularities and develops a smooth neck. 
As it is shown in [9]- [TO], the condition for bubbles to be at equal pressure implies that the integral 

where uj is the meromorphic extension of § dX to 1Z and (3 is a cycle connecting the bubbles, is a 
constant of the motion. Since at break-off (5 contracts to a point, it is obvious that a necessary 
condition for break-off is that II vanishes. 




The following pair of complex-valued functions were introduced in [10] to describe the bubble 
break-off near the breaking point 

X(X):=l(\ + S{\j), Y(\):=±(\-S(\j). (41) 



2i 

They analytically extends the Cartesian coordinates (X, Y) of the interface 7 

X = Re\, Y = ImX, A e 7, (42) 
and allow to write the evolution law (|40|) in the form 

d t Y(X) = -dxHX). (43) 



The analysis of [10] concludes that after the break-off the local structure of a small part of the 
interface containing the tips of the bubbles falls into universal classes characterized by two even 
integers (4n,2), n > 1, and a finite number 2n of real deformation parameters tk- By assuming 
symmetry of the curve with respect to the X-axis, the general solution for the curve and the 
potential in the (4n, 2) class are 

Y :=(-=%_=) y/(X-a)(Y-b), * = -y/(X - a)(Y -b), (44) 
V y/(X - a)(X -6)/© 

where a and b are the positions of the bubbles tips and 

In 

U(X,t):=J2tj+i xj+1 - (45) 
3=1 

Here the subscript © denotes the projection of X-series on the positive powers. Due to the physical 
assumptions of the problem, the function Y inherates two conditions for its expansion as X — > 00 

2n 00 v 

Y(x) = j2ti+i)tj + ix j +Y,j^- ( 46 ) 

3=1 j=0 

which determine the positions a, b of the tips. The conditions are 
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1. From (|38p <& — * — i log A as A — > oo. Hence (|43p implies that the constant term Yq in (|46p 
should be equal to t. 

2. The coefficient Y\ in front of X~ l turns to be equal to II, so that it must vanish for a break-off 

As it was shown in [10] , imposing these two conditions on (|46h leads to a pair of hodograph equations 
which arise in the dispersionless AKNS hierachy. However, from (j44H it is straightforward to see 
that these equations coincide with the hodograph equations (|24|) associated with the system of 
string equations (jlj) provided one sets 

X = z, Y = 2m-V z , $ = z -u-2p, 

(47) 

tj = 0, Vj>2n + 2; t = t l , x = — = 0. 

For instance, we observe that the evolution law (|43[) derives in a very natural form from the dl-Toda 
hierarchy. Indeed, from (|13|) and (|44|) we have 

p = - (z-u- 

so that (|4Tj) implies 

9 t y = 2d tl m(z) - 1 = 2d z (z)+ -l = 2d z (p + u)-l = -8 Z <S> = -d x *. 

In this way the integrable structure associated to the system of string equations (jl]) of the d2- 
Toda hierarchy manifests a duality between the planar limit of the Hermitian matrix model and 
the bubble break-off in Hele-Shaw cells. According to this relationship the density of eigenvalues 
p and the end-points a, b of its support in the Hermitian model are identified with the interface 
function Y and the positions of the bubbles tips, respectively, in the Hele-Show model. 
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